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1. Introduction 

The notion of stability in Dynamical Systems refers to dynamical behavior that persists 
under perturbation.^ By altering the nature of the persistence and the class of perturba- 
tions one obtains various forms of stability. These various forms of stability have proved to 
be extremely important throughout the history of dynamics. In perhaps the best known 
cases, KAM theory and structural stability, the dynamical stability involves small pertur- 
bations. However, dynamical persistence under large perturbations (in a restricted class) 
is often studied and has proved to be quite powerful. Large perturbation theories usually 
have a strong topological component. This is because behavior that persists under large 
perturbations must be very fundamental to the system, and the most fundamental aspect 
of a dynamical system is the topology of the underlying manifold. 

In applying stability results, one often begins with a model system whose dynamics are 
understood and then perturbs it. The stability theorems indicate which dynamics of the 
model system must be present in the perturbed system. This strategy often yields a great 
deal of information about the perturbed system that could not be otherwise obtained. In 
particular, it frequently provides a framework for the investigation the other dynamics 
present in the perturbed system. Since the dynamics of the model system must be present 
in all perturbed systems the model system may also be viewed as a dynamically minimal 
element in the allowed class. 

This paper surveys various results concerning stability for the dynamics of Lagrangian 
(or Hamiltonian) systems on compact manifolds. The main, positive results state, roughly, 
that if the configuration manifold carries a hyperbolic metric, i.e. a metric of constant 
negative curvature, then the dynamics of the geodesic flow persists in the Euler-Lagrange 
flows of a large class of time-periodic Lagrangian systems. This class contains all time- 
periodic mechanical systems on such manifolds. These results are given in Theorem 1 and 
Theorem 3 in Section 3. Complete proofs appear in [10]. Many of the results on Lagrangian 
systems also hold for twist maps on the cotangent bundle of hyperbolic manifolds. 

We also present a new stability result for autonomous Lagrangian systems on the 
two torus (see Theorem 5, Section 5) which shows, among other things, that there are 

^In other contexts "stability" may mean that smaU perturbations of initial conditions in (some region 
of) the domain of a fixed system give rise to small changes in asymptotic behavior. We will always use 
stability to refer to persistence of the global dynamics under perturbations of the entire system. 
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minimizers of all rotation directions. However, in contrast to the previously known [19] 
case of just a metric, the result allows the possibility of gaps in the speed spectrum of 
minimizers. Our negative result is an example of an autonomous mechanical Lagrangian 
system on the two-torus in which this gap actually occurs. The same system also gives us 
an example of a Euler-Lagrange minimizer which is not a Jacobi minimizer on its energy 
level. 

Our results generalize several theories that contain what may be viewed as stability 
results. The first is the ATibry-Mathcr theory. This theory shows that an area preserving 
monotone twist map of the annulus always has nicely behaved invariant sets with each 
rotation number. These invariant sets can be viewed as the remnants of the invariant 
circles of a minimal model, the rigid twist, or equivalently, the time one map of the full 
geodesic flow of the Euclidean metric on the circle. The Aubry-Mather theory is closely 
related to Hedlund's work on geodesies on the two torus ([19], cf [5]). Hedlund showed 
that for any Riemannian metric there are geodesies with all rotation directions, and thus 
the model system in this case is the Euclidean metric on the torus. In related work, closely 
connected to the hyperbolic manifold results here, Morse [32] showed that any metric on 
a higher genus surface has a collection of geodesies that "shadow" in the universal cover 
the geodesies of the hyperbolic metric. There are also generalizations to any dimensions 
and improvements of Morse's results due to Klingenberg [22], Gromov [17], and MacKay 
and Denvir [30]. Because we allow time dependent Lagrangians, these results do not imply 
ours. 

All these theories share the property that the orbits of the dynamical system under 
consideration correspond to extremals of a variational problem defined in the universal 
cover of the configuration space. The orbits that correspond to minima of the variational 
problem have special properties; they behave approximately like the solutions to the vari- 
ational problem associated with the model system. This enables one to take limits of 
minimizing orbit segments or minimizing periodic orbits in order to construct a large set 
of minimizing orbits on which the dynamics is similar to that of the unperturbed sys- 
tem. It is natural to study minimizers in the perturbed systems because all orbits of the 
unperturbed system are minimizers. 

There is a simple heuristic connection between these different theories. In the Aubry- 
Mather theory, Aubry's Fundamental Lemma [4], [29], states that minimizers for a twist 
map are ordered like orbits under a rigid rotation ( i.e. like orbits of the time-1 map of the 
Euclidean geodesic flow on the circle). This is easily seen to imply that such orbits have 
a rotation number and the rotation number of a limit of such orbits is the limit of the 
rotation numbers. If {xk}kez is such an orbit and it has a rotation number uj = limxk/k, 
then one checks that \xk — xi — {k — 1)lo\ < 1 for all k,l & Z. This immediately implies 
that 

\uj\\k -l\-l<\xk- xi\ < \uj\\k + 

This, in Gromov's language, says that the well ordered orbits are quasi- geodesies. It turns 
out that most of the old or more recent Riemannian geometry results on the stability of 
the hyperbolic geodesic flow can be proved using the fact that minimal geodesies for any 
metric are quasi-geodesics for the hyperbolic metric. In our work on hyperbolic manifolds 
[10] that we survey here, we also use a limiting argument, and our central Proposition 2 
(which is rather trivial in the autonomous case) states that Euler-Lagrange minimizing 
segments of a given average speed are quasi-geodesics. The proof of this proposition uses 
techniques that go back to Aubry's proof of his Fundamental Lemma, which have a parallel 
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in Riemannian geometry under the guise of curve shortening arguments. We should remark 
here that the property of being a quasi-geodesic is a rather weak regularity property and 
even though it is satisfied by all minimizers for our large class of systems on any manifolds, 
it is not sufficiently strong to make the limiting argument work on manifolds that do not 
support a hyperbolic metric (except in cases of very low dimension) . We illustrate this in 
Section 3.3 with the Hedlund example on the three torus. 

We were greatly influenced by and used many techniques of the recent work of Mather 
[27], [28] (see also [25]), which attempts, among other things, to generalize the Aubry- 
Mather theory to higher dimensions. Moser showed that convex, time-periodic Lagrangian 
systems on compact manifolds are a natural generalization of twist maps in that twist 
maps arc always the time-1 maps of such Lagrangians on the tangent space of the circle 
([33]). Denzler [12] followed through this philosophy and gave a proof of the Aubry-Mather 
Theorem in the larger context of time-periodic, convex Lagrangians on the circle. 

The class of Lagrangians we consider is almost the same as Mather's. One fundamental 
difference, however, is that Mather uses minimizers not in the universal cover as is done 
here in the hyperbolic case, but rather in the universal free Abelian cover (the cover 
with deck group H{M;Z) /torsion). If the fundamental group of the configuration space is 
torsion- free Abelian (eg. a torus) this cover is the universal cover, but in general, it is much 
smaller than the universal cover. Whereas Mather's theory works on any compact manifold, 
in the special case of hyperbolic manifolds, much information is lost. Indeed, many orbits 
of the hyperbolic geodesic flow are not minimizers in the sense of Mather (see Example 5). 
Hence there is no chance in finding corresponding orbits in a perturbed system by looking 
for such minimizers. On the other hand, in our results on the torus presented here, we 
use the full strength of Mather's results (which we review for the reader's convenience), 
in particular the Lipschitz graph property of his generalized Mather sets. 

It is important to remark, especially in the context of this conference, that Mather's 
stated goal in his recent work is not so much to prove stability results for their own sake but 
to use his Mather sets (the remnants he finds of the dynamics of the unperturbed system) 
as a stairway to (generic) diffusion. He accomplished this program in [26] for twist maps 
(which are time-1 maps of time periodic Lagrangians on S^), and gets partial results for 
the general case in [28]. We hope that our work might help to reach this important goal. 

This paper is organized as follows. Section 2 introduces the Lagrangian setting and 
several fundamental lemmas are given. In Section 3, we present and outline the sketch 
of our two theorems on time periodic Lagrangians on hyperbolic manifolds. The first 
(Theorem 1 in Section 3) states that you can shadow any geodesic of the hyperbolic 
metric by a minimizer of the Lagrangian at uniform bounded distance. This theorem also 
holds for symplectic twist maps in T*M. The second (Theorem 3) concerns large invariant 
sets for the Euler Lagrange flow which are built from these shadowing minimizers, and 
shows that the dynamics on the sets is semiconjugate to the hyperbolic geodesic flow. At 
the end of the section we remark on why this scheme of proof does not work on the three 
torus. 

In Section 4, we briefly review Mather's recent work on minimal measures, and illustrate 

it by examples. In Section 5, we apply Mather's theory to autonomous systems on the two 
torus, and give a fairly complete description of the rotation set of minimizers in this case 
(Theorem 5, Section 5). In particular, we find orbits of all rotation directions, and infinitely 
many of different average speeds in each direction. 

In Section 6, we show that the gaps in speed that are allowed in Theorem 5 actually 
occur in natural, mechanical systems. 
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2. Preliminaries 

In this section we introduce notation and recall some basic results needed in the sequel. 
For a thorough discussion of Lagrangian systems and minimizers the reader is urged to 
consult Mather [27] and Mane [25] . We also indicate how to translate the different notions 
to the setting of symplectic twist maps. 



2.1. LAGRANGIAN SYSTEMS. 

The main objects in the Lagrangian fornnilation of mechanics are a configuration manifold 
M and a real valued function called a Lagrangian defined on the tangent bundle TM. The 
configuration spaces of interest here are closed manifolds M with a fixed Riemannian 
metric g. The induced norm on the tangent bundle is denoted ||f ||. We consider time- 
periodic systems determined by a C^-Lagrangian L : TM x §^ — > R. The basic variational 
problem is to find curves 7 : [a, 6] — M that are extremal for the action 

^(7)= [\{-f,^,t)dt 

J a 

among all absolutely continuous curves /? : [a,b] — M that have the same endpoints 

/3(a) = 7(0), /3(6)=7(&). 

Under appropriate hypothesis {eg. 7 is C^), such a 7 satisfies the Euler-Lagrange second 
order differential equations 

If (7(*),7W,0-f (7(t),7(t),t) = 0. 

Using local coordinates these equations yield a first order time-periodic differential equa- 
tion on TM, and thus in the standard way, a vector field on TM x S^. Since TM x is 
not compact it is possible that trajectories of this vector field are not defined for all time 
in R and thus do not fit together to give a global flow {i.e. an M-action). When the flow 
does exist, it is called the Euler-Lagrange or E-L flow . 

We will require that Lagrangians satisfy certain hypotheses. 

Hypothesis: 

L is a C'^ function L : TM xS^ that satisfies: 
(a) Convexity: is positive definite. 

(h) Completeness: The Euler-Lagrange flow determined by L exists, 
(c) Superquadratic: There exists a C > so that L{x,v,t) > C||f p. 
or 

(c') Superlinear: — ^ 00 when \\v\\ — -|-oo. 

We will refer to (a), (6) and (c) as Our Hypotheses and (a), (6) and (c') as Mather's 
Hypothesis. Note that our hypotheses are a little stronger than Mather's (adding a constant 
to the Lagrangian doesn't change the E-L flow). 

Example 1 : Mechanical Lagrangians. As pointed out by Mane, Mather's Hypothesis 
( and hence ours ) are satisfied for mechanical Lagrangians, i. e. those of the form 

L{x,v,t) = ^\\vf - V{x,t), 

where the norm is taken with respect to any Riemannian metric on the manifold. (In fact, 
one may allow the norm to vary with time, under some conditions, see [25], page 44)- 
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2.2. MINIMIZERS 

Of particular interest here are extremals of the variational problem that minimize in the 
following sense. If M is a cover of M, L lifts to a real valued function (also called L) defined 
on TM X S^. A curve segment 7 : [a, 6] — M is called a M-minimizing segment or an M- 
minimizer if it minimizes the action among all absolutely continuous curves /? : [a, 6] — M 
which have the same endpoints. 

A fundamental theorem of Tonelli implies that if L satisfies Mather's Hypotheses, 
then given a < b and two distinct points Xa,Xb G M there is always a minimizer 7 with 
7(a) = Xa and 7(6) = Xf). Moreover such a 7 is automatically and satisfies the Euler- 
Lagrange equations (this uses the completeness of the E-L flow). Hence its differential 
dl{t) = (7(i), 7(t)) yields a solution {d'y{t),t) of the E-L flow. A curve 7 : M — M is called 
a minimizer if 7|j^ is a minimizer for all [a,b] C M. When the domain of definition of a 
curve is not explicitly given it is assumed to be M. 

As noted in the introduction, Mather [27] and Marie [25] use M-minimizers where M 
is the universal free Abelian cover. The universal cover (which we denote M from now on) 
is used here. If 7 is an M-minimizer, we will simply say it is a minimizer. 

Our main task is to get control of the speed and geometry of M-minimizers. Given a 
smooth curve 7 : [c, d] — M and a segment [a, b] C [c, d], the average displacement in the 
cover over the interval [a, b] is measured by 

b — a 

where d is the topological metric on M constructed from the lift of the given Riemannian 
metric g. 

The fact that L is superlinear or superquadratic implies some very useful simple esti- 
mates on the average action of minimizers. These estimates are essentially in Mane [25] 
and Mather [27]. 

Lemma 1 Given a Lagrangian L satisfying Mather's Hypothesis, there are functions K 
C^^^,K I— ^ C^"^ both increasing to infinity such that if ^ is a minimizer and S{'y;a,b) = 
K, then 

CT"K < / L(7,7,t)dt < C^«^K. 

b- aJa 

In particular, if L = ^Hxjl^ — V{x,t), we have: 

1 

Vmax < 7 / L{l,'y,t)dt < — Vmin- 

Z a J a z 



Proof Referring the reader to [25] or [10] for a proof of the first, general statement of the 
lemma, we give a proof of the mechanical case. We first estimate the lower bound: 



h 

L{-f,'y,t)dt 



b - aJa 

1 1 f y 
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- ^-V 

— 2 max • 

The inequality on the second line follows from Cauchy-Schwarz: J J > {J fg)"^, setting 
/ = ||x|| and (7 = 1. 

To get the upper bound, let T : [a, 5] — > M with 7(a) = T(a) and 7(6) = T{b) be a length 
minimizing geodesic segment with respect to the given metric. Then ||r|| = p(r;a,b) = 
5(7; a, b) = K. Thus, since 7 is a minimizer, 

^(7|[„,,]) < ^(r) < j^^-Vmrn<[—- Vmin] {b - a), 

yielding the upper bound. □ 



2.3. EXACT SYMPLECTIC TWIST MAPS 

For more details on symplectic twist maps, the reader is referred to [16] or [23] (see also 

[7] and [21]). An exact symplectic twist map F is a map from a subset U of the cotangent 
bundle of a manifold N (which we allow to be noncompact) into U, which comes equipped 
with a generating function S : N x N ^M. that satisfies 

F*{p dx)-pdx = P dX -pdx = dS{x, X), (1.2) 

where {X,P) are the coordinates of F{x,p) (this can also be written in a coordinate free 
manner) . 

Because the one-form P dX — p dx in (1.2) is exact, one says that F is exact. Note 
that taking the exterior differential of (1.2) yields dP A dX = dp A dx, and so any exact 
F is also symplectic, i.e. it preserves the standard symplectic form. The fact that S is 
expressed using the coordinates {x,X) instead of {x,p) is the twist condition. Given S, 
one can retrieve the map (at least implicitly) from p = —^andP = This can be done 
globally (i.e. U = T*N) only when N is diffeomorphic to a fiber of T*N, for example when 

is the covering space of the n-torus or of a manifold of constant negative curvature. 

The variational problem for Lagrangian systems translates into a discrete variational 
problem for twist maps: the role of curves in the continuous setting is taken by sequences 
of points ("integer time curves"), and the action of a finite sequence x = {xn, ■ ■ ■ ,Xm} is 
given by W{x) = J2^~^ S{xk,Xk+i). This corresponds closely to the continuous setting 
when the exact symplectic twist map F is the time-one map of an E-L flow. In this case, 
S{x, X) = L{x, X, t)dt, where x{t) is the minimizer over the interval [0, 1] with endpoints 
X and X. 

In direct correspondence to Lagrangian systems, critical points of W (with fixed time 
and configuration endpoints) correspond to orbits of F (this is closely related to the 
method of broken geodesies in Riemannian geometry). Action minimizers are sequences 
that minimize W over any of their subsegments. The natural growth condition on the 
generating function 

S{x,X) >Cdist^ix,X), 

implies the analog of Tonelli's theorem: minimizers always exist between any two points 
over any given (integer) interval of time. Moreover, there is an exact analog of Lemma 
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1.1: the average action of minimizers is bounded below and above by functions of the 
average displacement. The proof is virtually identical to the continuous time case, replacing 
geodesies with orbits of the time-one map of the geodesic flow. 

Example 2 Let M be T" or a closed hyperbolic manifold, and let N = M be the universal 
cover R" or H", respectively. On the covering space, define the generalized standard map 
using its generating function M x M — > R, 

S{x,X) = ^disf{x,X) + V{x) 

where the distance dist is induced by the Euclidean metric in R", or the hyperbolic metric 
on IP, and V{x) is 7ri(M)-equivariant, i.e. it descends to a function on M. A short 
argument shows that one can use the relation (1.2) to solve for {X,P) in terms of {x,p) 
and thus obtain an exact symplectic twist map on T*M that, in turn, induces a map on 
T*M (also called a twist map). For more general examples, cf [16]. 

In certain cases the twist map theory overlaps with the continuous theory. If a twist 
map / of T*T" has a generating function that is super quadratic in \\X — x\\, the mixed 
partial 8128 is symmetric, and for some a > satisfies the convexity condition 

< di2S{x,X).v,v > < -a\\v\f 

uniformly in (x,X), then F is the timc-onc map of an E-L flow derived from a one- 
periodic Lagrangian that is superquadratic in the velocity. Moser [33] gives the proof in 
the case n = 1. Bialy and Polterovitch remark in [11] that Moser 's proof goes through 
in the case n > 1. This is not quite so, but they subsequently obtained a different proof 
(personal communication) . Note that the generating function for the generalized standard 
map satisfies these hypothesis. 



2.4. JACOBI MINIMIZERS VS. E-L MINIMIZERS 

In this section, we consider autonomous mechanical Lagrangian of the form L{x,x) = 
(l/2)||xf -y(x) where < and || comes from a Ricmannian metric on the manifold 
M. In this case there is a geometric way to look for E-L minimizers on a given energy 
level. Recall that, in the tangent bundle coordinates, the energy (Hamiltonian) is given by 
H{x,x) = L{x,x) + 2V{x) = {l/2)\\xf + V{x) 

The norm on the tangent bundle coming from the Jacobi metric with energy E > V„iax 
is \JE — 1/(x)||i;||. The important fact is that geodesies of the Jacobi metric considered as 
curves in M are always the projection to M of some solution curve with energy E of the 
Hamiltonian (or E-L) flow, and conversely, all such projections are geodesies of the Jacobi 
metric. However, the parameterizations of the geodesies and the solutions will usually be 
different, (eg see Arnol'd [3] or Abraham and Marsden [1]) 

Put in other language, this means that the extrema of the corresponding integrals 
coincide in some sense. What is of importance here is whether the minimizers coincide. 
Let us define the notion of minimizer more carefully. 

Definition 1 1. A curve 7 : [a,b] ^ M is called an E-L minimizer if for all absolutely 
continuous /3 : [o, 6] — > M with (3{a) = 7(a) and P{b) = 7(6), 

t L{^,^)dt< f\{l3,^))dt 

J a J a 
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2. A curve 7 : [a,b] ^ M is called a Jacobi minimizer for energy E if for all absolutely 
continuous {3 : [ap, bp\ — > M with P{ap) = 7(a) and P{bp) = 7(6), 



Note that being an E-L minimizer is a property of both the path and the parameter- 
ization. On the other hand, being a Jacobi minimizer just depends on the path since the 
integral in its definition is 1-homogeneous in the derivative. 

The fact that solutions to the E-L equations are extrema of the integral in (1) is usually 
called Hamilton's Principle. On the other hand, The Principal of Least Action states that 
if 7 is a solution to the E-L equations in M, then among the curves (3 with /9(a^) = 7(a) 
and P{b/s) = 7(6) which satisfy the constraint H{l3,$) = E, the curve 7 will extremize 



the integral / . Now if H{(3,f3) = E, then = y'2{E - V{fj))\\p\\, and so the 



Principle of Least Action is the same as the variational formulation for geodesies with 
the Jacobi metric as in (2). (The reader is warned that there is a great deal of variance 
in the literature regarding the meaning of the phrases "Principle of Least Action" and 
"Hamilton's Principle".) 

Proposition 1 If j : [a,b] ^ M is a Jacobi minimizer with energy E > Vmax, then there 
exists a re-parameterization 71 with i/(7i,7i) = E that is an E-L minimizer. 

Proof Begin by noting that 



and equality occurs if = E. Indeed, setting a = y/E-V,b= \\$\\/V2, the above 

inequality is just -\- b^ > 2ab, and a = b means H = E. 

Assume that 7 is a Jacobi minimizer on the energy level {H = E}, and so we may 
parameterize 7 so it is a solution of the E-L equations. Let this parameterization be 
71 : [a, b] M. Now let /? be a test path as in the definition of E-L minimizer. Using the 
inequality of the previous paragraph and the fact that 7 is a Jacobi minimizer we have 



and so 71 is also a E-L minimizer. □ 

Remark 1 Somewhat surprisingly, the converse of the proposition is false; there can be 
E-L minimizers which are not Jacobi minimizers on their energy levels. What basically 
happens is that there are two E-L minimizers connecting the same two points in configu- 
ration space, but they minimize over different time intervals. They have the same energy 
but different Jacobi lengths, and thus they both cannot be Jacobi minimizers. An example 
is given at the end of Section 6. 

As noted above. The Principle of Least Action is the same as the variational formulation 
of Jacobi geodesies. Thus there are minimizers of / L with fixed time- and space-endpoints 




LiP,$) + E>^2{E-Vmm 




<j^vHE-v{mm 
<j^m$)+E 
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that are not minimizers of the action integral with fixed space-endpoint and constant 
energy constraints. 

Since a solution of the E-L equations is an cxtrcma and not necessarily a minimum 
of the action integral, many authors have adapted the terminology "Principle of Extreme 
Action". Our example gives another argument for this terminology: while the extrema of 
the Euler Lagrange and Action problems coincide, their minimizers may not. 

3. Stability Results for Lagrangian Systems on Hyperbolic Manifolds 

We consider here a closed hyperbolic manifold of arbitrary dimension n, equipped with a 
metric of constant negative curvature. All such manifolds M have the Poincare n-disk H" 
as universal covering: M = H", with the canonical hyperbolic metric. Geodesies for this 
metric in are arcs of (Euclidean) circles perpendicular to the sphere at infinity (the 
Euclidean unit sphere). The complete proofs of Theorems 1 and 3 appear in [10]. 

3.1. SHADOWING GEODESICS WITH E-L MINIMIZERS 

Our first result gives the first way of formalizing the notion that Lagrangian systems satis- 
fying our hypothesis (see Section 2.1) on hyperbolic manifolds are at least as complicated 
as the geodesies of a hyperbolic metric. Given a hyperbolic geodesic in the Poincare Disk 
H", the theorem asserts that there are minimizers of the Lagrangian system that are a 
bounded distance away and have a variety of approximate speeds. Recall that 5{j;a,b) 
means the average displacement in H" over the time interval [a,b], i.e. the distance from 
7(a) to 7(6) divided hy b — a (see Section 2.1). 

Theorem 1 Let (M, g) be a closed hyperbolic manifold. Given a Lagrangian L which sat- 
isfies our hypothesis there are sequences ki,Ki,Ti in depending only on L, with ki 

increasing to infinity, such that, for any hyperbolic geodesic T C H" = M , there are min- 
imizers 7j : M — > M with dist{'yi,r) < Ki, 7i(±oo) = r(±oo), and ki < 6{'y;c,d) < Aij+i 
whenever d — c>Ti. 

Remark 2 This theorem also holds for twist maps of T*M provided their generating 
function has super quadratic growth, eg. the generalized standard map with generating 
function S{x,X) = {l/2)dist^{x,X) -\- V{x). The proof is identical to that of the La- 
grangian case, using the dictionary given in Section 2.3. 

Proof (Sketch) Fix an oriented geodesic with a given parameterization by arclength 

r : M ^ H" and a. K > Kq with Kq as in Proposition 2 below. Let -/n ■ [-N, N] 
be a minimizing segment with 7Ar(— iV) = r{—KN) and 7jv(-^) = ^{KN), and thus 
S{lN', —N, N) = K. We would like to be able to take a limit of the curves jn's and prove 
that the limit is a E-L minimizer with the same endpoints at infinity as P. In order to do 
this, we use Gromov's theory of quasi-geodesics. 

Definition 2 Given A > 1 and e >0, a curve 7 : R — ^ H" or a curve segment 7 : [a, b] — ^ 
H" is called a {X,€) -quasi-geodesic if 

X-\d - c) - e < d{j{c), 7(d)) < A(d - c) + e 

for all [c,d] in the domain of^. 
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The next theorem, often called "Stability of quasi-geodesics" , gives the most important 
property of quasi-geodesics. It is true in the broader context of what are usually called 
(5-hyperbolic spaces, but we just state the result in the context needed here. Given two 
closed subsets X,Y C H", d(X,Y) denotes their Hausdorff distance as induced by the 
hyperbolic metric. For a proof and more information see [15]: 

Theorem 2 Given A > 1 and e > 0, there exists a n > so that whenever j is a (A,e)- 
quasi- geodesic segment in and Tq is the geodesic segment connecting the endpoints of ^, 
then di^jFo) < k. If ^ is a {X,€) -quasi-geodesic, then j{oo) and ^{—oo) exist and further, 
ifT is the geodesic connecting 7(00) and 7(— 00), then d{'y,T) < k. 

Proposition 2 There exists Kq > depending only on L and A > 1 and e > depending 
only on K and L so that whenever 7 : [a, 6] — > M «s a minimizing segment with (5(7; a, 6) = 
K > Kq and b — a EN, then j is a (A, e)-quasi-geodesic segment. 

Thus, our proposition implies that the 7iv defined in the first paragraph of the proof of 
Theorem 1 are (A, e)-quasi-gcodesics for the same (A,e). Theorem 2 implies that the 7Ar's 
stay at uniformly bounded distance from T. This enables us to take a (pointwise) limit of 
the 7iv's and get an E-L minimizer 7 which has same endpoints at infinity as T and stays 
at bounded distance from F. 

The proof of Proposition 2 is somewhat technical. We just sketch the essential features. 
We have that our minimizer 7 satisfies (5(7; a, b) = K. Since 5{^\ a, b) = ^^iM^iilM) ^ this 
translates to dist{'y{a),'y{b)) = K{b — a). From this estimate at the endpoints of 7, we 
want to derive similar estimates for any interval [c, d] C [a, 6]. Ideally, we would like to 
prove that 

K' < (5(7; c, d) < K" (1) 

for K' , K" only depending on K. This would immediately prove Proposition 2. The second 
inequality follows from a results of Mather: (5(7; a, 6) = K ^ \\^\\ < K" where K" only 
depends on K. The first inequality, proved in [10], is only true when (d — c) > N, for 
some N only depending on K. But this restriction only adds an e in the quasi-geodesic 
estimate. Note that, in the autonomous case, these estimates are trivial. 

To prove the second inequality in Formula 1, Mather uses a surgery argument of a 
type which is widespread in the Aubry-Mather theory. If 7 was going very fast on a 
subinterval [c,d], then one could cut 7|[^^j out and replace it by a minimizer between the 
same endpoints but on a greater interval of time, thus reducing the average speed on this 
interval. To make up for lost time one does the same kind of surgery on another judiciously 
chosen interval on the curve. The action of this new curve 7* : [a, 6] ^ M is then estimated 
using the average action vs. average speed estimate of Lemma 1 of section 2.2, and shown 
to be less than that of 7, which is absurd. Our proof of the first inequality uses similar 
techniques. 

To finish the sketch of the proof of Theorem 1, we note that, in our construction 
of the minimizers 7jv's, we can vary K. Because of the bounds on the average speed 
on subintervals we explained above, this implies that, if Ki is chosen to be increasing 
sufficiently fast to 00 with i, the corresponding limiting 7^ are distinct and of average 
speed increasing to infinity. 
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3.2. SEMICONJUGACY WITH THE GEODESIC FLOW 

One way to formulate the fact that the perturbed system is at least as complicated as 
the model system (i.e. the dynamics of the model systems don't go away) is to show 
that the perturbed system always has a invariant set that carries the dynamics of the 
minimal model. More precisely, one shows that there is a compact invariant set that is 
semiconjugate to the minimal model (this strategy is common in more topological theories, 
see [8]). MacKay and Denvir [30] have recently extended Morse's results to the case with 
boundary and proved a result giving this semiconjugacy. Also Gromov [17] and others 
have done this in the case of geodesic flows. 

Theorem 3 Let M be a closed hyperbolic manifold with a hyperbolic metric g with geodesic 
flow gt- Given a Lagrangian L which satisfies the Hypotheses of Section 2 with E-L flow 4>t, 
there exists a sequences ki and with ki increasing to infinity, and a family of compact, 
(f)t-invariant sets Xi C TM x S-^ so that for all i, {Xi,(j)t) is semiconjugate to {TiM,gt) 
and ki < 6{(j)t{x); 0, T) < ki+i, whenever T > Ti and x ^ Xi. 

Note that the geodesic flow of a hyperbolic metric is transitive Anosov and thus is 
Bernoulli, has positive entropy, etc. Thus Theorem B implies that the E-L flow is always 
dynamically very complicated. 

We first state precisely what we mean by semiconjugacy : 

Definition 3 Two flows {X,(/)t) and (Yjipt) cl^c- said to be semiconjugate (or sometimes 
orbit semi- equivalent) if there is a continuous surjection f : X ^ Y that takes orbits 
of (pt to those of ipt preserving the direction of the flow, but not necessarily the time 
parameterization. Note that f is a local homeomorphism when restricted to an orbit of (f)t, 
but f may take many orbits of (pt to the same orbit ofipt- 

Given a K and a geodesic F, we have constructed a (A, e)-quasi-geodesic 7^ which 
shadows T. The differential of such a curve is in TM. Project all these differential curves 
down to TM and take the closure of this set in TAI. This gives a compact invariant set Qk 
for the E-L flow because the velocities of the uniformly bounded. Finally, define 

Qk to be the set of all possible lifts to TM of all the points in Qk ■ 

We will show that the E-L flow restricted to Qk is semiconjugate to the geodesic flow. 
Qk is a set made of (A, e)-quasi-geodesics, for some fixed (A, e). Each such quasi-geodesic 
shadows a unique geodesic, by Theorem 2, and each geodesic is shadowed by at least one 
quasi-geodesic in Qk, by Theorem 1. Hence we have a well defined application z i-^ 
which, to a point z in Qk, makes correspond the unique geodesic C H" that the 
quasi-geodesic to which z belongs shadows. 

We then project z on dV^ in the following fashion. Fix a parameterization by arclength 
for each geodesic F. Project 7r{z) G on F^ via the orthogonal projection (i.e. by drawing 
the unique geodesic through tt{z) which is perpendicular to F^) and get a point Tz{s{z)) 
on F^. Define 

a(z) = (F,(s(z)),f,(s(z)). 

Theorem 1 implies that a : Qk TiM is onto. It is not too hard to see that a is also 
continuous, and equivariant, i.e. it descends to a continuous map Qk TiM. Unfortu- 
nately, it is not necessarily injective when restricted to an orbit of the E-L flow in Qk- 
This is remedied using an averaging technique due to Fuller [14]. 

Given z = z(0) G Qk and its orbit {z[t),t G R} under the E-L flow, let a{z,t) = 
s{z{t)) — s{z{0)), where, s{z) is as above. If a{z,t) were positive for all t, we would be 
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done: a would be injective along z{t). Because we cannot assume that, we use the fact that 
a{z,t) is positive for large t and "average" a over the interval [0,t]. Given a > define 

(Toiiz) = a(z) -\ — / a(z,t) dt. 
a Jo 

Informally, aa{z) is the average value of a over the orbit segment 2;([0, a]). 

Now since for every z € Qk we have that uj{z) = r2;(oo), it follows that for each z there 
is an so that a{z,az) > 0. Since Qk is compact, we may find an a with a{z,a) > 
for all z G Qk- Let a = aa- As before, we write a{z) = {Tz{s{z)),Tz{s(z)). Now a is 
clearly continuous, equivariant, onto and takes orbits to orbits. We then show [10] that it 
is injective on orbits z{t) in Qk by showing that for any t > 0, s{z{t)) — s(z(0)) > 0. This is 
done by using the fact that a is an additive cocycle, i.e. a{z, ti+t2) = a{z, ti)+a{(pti {z),t2), 
for all ti,t2. By repeating this whole process by setting K = Ki, where Ki is as in Theorem 

1, we get distinct Qk^ which arc semiconjugate to the geodesic flow, and such that the 
average speed of a E-L orbit on Qxi goes to oo with i. This finishes the sketch of our proof 
of Theorem 3. 

3.3. WHY THE LIMIT ARGUMENT DOES NOT WORK IN 

Quasi-geodesics can obviously be defined on any manifold (in fact, on any metric space). 
Proposition 2 is valid on any compact manifold, not just those that support a hyperbolic 
metric. The reason why the proof above does not apply to M = must therefore be that 
quasi-geodesics on do not satisfy Theorem 2. Thus a sequence of (A, e)-quasi-geodesic 
segments with endpoints on a fixed geodesic F do not have to stay at bounded distance of 

r. 

As a very simple example, one can look at M? with the Euclidean metric and take the 
geodesic F to be the line y = x. Take the sequence of "corner" curves going between the 
points {—N, —N) and (N, N) of F by first following the line x = —N upward at unit speed 
until it reaches y = N, and then follow that line to (N,N). It is easy to check that the 
77v's are quasi-geodesics segment for a (A, e) uniform in TV. But obviously d{-~fi\f,T) — oo. 
It is not too hard either to sec that one can adapt this argument in Hedlund's example 
(sec Section 4.2) in to show how our limiting scheme would fail in finding a minimizer 
of rotation vector (1,1,0) there: in this case, the minimizers 7Ar's between the points 
{—N, —N, 0) and {N, N, 0) of the (Euclidean) geodesic F given by the line y = a;, ^; = of 

would be corner curves that follow the tubes, jumping tubes at the beginning, corner 
and end of the curve. Again the -fjsf are quasi-geodesic, with uniform (A, e) by Proposition 

2, and again their distance to F goes to infinity. Note that 'Jn ^ Fa; U Fj^ in the Hausdorff 
topology, where F-^, Ty are the periodic orbits projection of the tubes parallel to the x and 
y axes. This is the Mather set for the vector of rotation (1, 1, 0) (sec next section). 

A very similar phenomenon happens in Example 5 of Section 4.2. In the case considered 
there the jn with "corners" show up in the universal free Abelian cover M. 

4. A Quick Review of Mather's Theory of Minimal Measures 
4.1. THEORY 

For a more detailed exposition the reader is urged to consult Mather [27] or Mane [25]. 
Mather also gives a very nice survey of this theory in the beginning of [28] . Given a E-L 
invariant probability measure with compact support fi on TM x one can define its 
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rotation vector p{iJ,) as follows: let Pi, P2, ■ ■ ■ , Pn he a basis of H^(M) and let Ai, . . . , A„ be 
closed one-forms with [Aj] = f5i in DeRham cohomology.^ The reader uncomfortable with 
homology may read through this section thinking of the case of M = T", with angular 
coordinates {xi, . . . ,x„), and taking [Aj] = [dxi], as a basis for H^{T"') ~ Hi{T^) ~ R"'. 
Define the i*^ component of the rotation vector as 



Note that this integral makes sense when one looks at A^ as inducing a function from 
TM X S-*^ to M by first projecting TM x onto TM, and then treating the form as a 
function on TM that is linear on fibers. The rotation vector does depend on the choice 
of basis /3j, but because the one forms are closed, pi{p) does not depend on the choice 
of representative Aj with [Aj] = Since the rotation vector is dual to forms, it can be 
viewed as an element of Hi{M,M). In the case M = T", one can check that, for 7 a generic 
point of an ergodic measure fj,, the usual rotation vector of 7 coincides with that of fj,: 

Pi{-f) = lim ^^^^ = lim / dxi = [ dxidp = pi{p) 

6— a— »oo — a 6— a— »oo — O J dri\^a b] 

where second equality uses the Ergodic Theorem. (Again, dxi is seen as a function TM x 
^ M). If M is a general compact manifold, one can define the rotation vector of a curve 
7 : M ^ M by pj(7) = limf,_a-+oo Sd-y\y^ ^ ^^e limit exists. As before, if 7(0) is a 

generic point for an ergodic measure ji, this rotation vector does exist and coincides with 
that of p. 

Next we define the average action of a E-L invariant probability on TM x S-*^ by 



A{p) = j Ldp, 



which, when p is ergodic, we can relate to the average action along p-a.e. orbit 7 by 

rb 



A{p)= lim / L{j,j)dt. 

0— a— »oo — a Ja 



The set of invariant probability measures, denoted A^l is a convex set in the vector 
space of all measures^ and the extreme points of A^l are the ergodic measures [24] . Now 
consider the map Ml —>■ Hi{M) x R given by: 

/^i-^ (p(/^),^(/^)) • 

This map is trivially linear and hence maps Ml to a convex set Ul whose extreme points 
are images of extreme points of A^l, i.e. images of ergodic measures. Mather shows, by 
taking limits of measures supported on long minimizers representing rational homology 
classes, that for each uj, there is p such that p{p) = uj and A{p) < cx).^ Since L is 

^When homology and cohomology coefficients are unspecified they are assumed to be M, so the notation 
Hi{M) means Hi{M;W), etc. 

^It is also compact for the weaJi topology if, as Mather does, one compactifies TM. 

''The impatient reader may be tempted to proclaim, from this fact, the existence of orbits of all rotation 
vectors. Alas, one can only deduce the rotation vector of orbits from that of a measure when the measure 
is ergodic... 
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bounded below, the action coordinate is bounded below on Ul- Hence we can define a 
map /? : Hi{M) ^ M by 

f]{uj) = inf{^(^) I fi G Ml, p{l^) = c^}, 

which is bounded below and convex; the graph of /? is the boundary of Ul- 

We say that a probability measure /x G TWl is a minimal measure if the point {p{lJi)i A{ii)) 
is on the graph of (5. Hence, an extreme point (cj, of graph{f3) corresponds to at least 

one minimal ergodic measure of rotation vector to. It turns out that if p is minimal, p-a.e. 
orbit lifts to a E-L minimizer in the covering M of M whose deck transformation group 
is Hi{M;Z) /torsion {i.e. the universal cover when M = T"). Conversely, if /x is an er- 
godic probability measure whose support consists of M-minimizers, then /x is a minimal 
measure. 

Hence, each time we prove the existence of an extreme point {uj,P{lo)), we find at least 
one recurrent orbit of rotation vector u which is a M -minimizer. 

Another important property of (3 is that it is superlinear, i.e oo when ||a:|| 

GO. Since we will need the estimate later, we motivate this in the simple case where 
= lll^lP ~" V{x) and || • || comes from the Euclidean metric on the torus. If /x is any 
invariant probability measure, then 

|2 



A{n) = j Ldll> j Vraaxj djl 

^max (2) 



1 

> - 



2 

X dfi 



2 

= \\P{IJ)\^ -Vmax 

where we used the Cauchy-Schwarz inequality for the second inequality. 

The super linearity of (3 implies the existence of many extreme points for graph{P) 
(although in most cases still too few, as we will see at the end of this discussion). Indeed, 
/3's superlinear growth implies that its graph cannot have flat, or linear domains going to 
infinity. Any point (a;,/3(a;)) is part of at least one linear domain of graph{/3), which we 
call Sc (we suppress the dependence of conuj). Here, the index c denotes the "slope" of the 
supporting hyperplane whose intersection with Ul is exactly the convex and flat domain Sc 
(c can be seen as element of first cohomology). Let Xc be the projection on Hi{M) of Sc- 
Then the X^s are compact and convex domains which "tile" the space Hi{M). Extreme 
points of Xc are projections of extreme points of Sc- Hence there are infinitely many such 
extreme points, and infinitely many outside any compact set. Their convex hull is Hi(M), 
and in particular, they must span Hi{M) as a vector space. Since these extreme points are 
the rotation vectors of minimal ergodic measures, wc have found that there always exist 
at least countably many minimal ergodic measures and at least n = dim Hi (M) of them 
with distinct rotation directions. We will see in Hedlund's example that this lower bound 
can be attained. 

Finally, the generalized Mather sets arc defined to be Mc = support{Aic), where Aic 
is the set of minimal measures whose rotation vector lies in Xc- Let n : TM x §^ M xS^ 
denote the projection. Mather's main result in [27] is the following theorem. 

Theorem 4 (Mather's Lipschitz Graph Theorem) For all c G H^{M), Mc is a 
compact, non-empty subset of TM x The restriction of tt to Mc is injective. The 
inverse mapping Tr~^ : 7r(Mc) Mc is Lipschitz. 
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In the case M = T", Mather proves that, when they exist, KAM tori coincide with the 
sets Mc, and that they are in the closure of these sets (see also [21] for some related results). 
In a sense, Mather's theorem generalizes Birkhoff 's theorem on invariant curves of twist 
maps, which says that such curves have to be graphs (see [20], [11] for more straightforward 
generalizations of this theorem). The proof of the Lipschitz Graph Theorem (see [27] or 
[25]), which is quite involved, uses a curve shortening argument: if curves in 7r(Mc) were 
too close to crossing transversally, one could "cut corners" and, because of recurrence, 
construct a closed curve with lesser action than Amin- 

Remark 3 An important special case is that of autonomous systems. In this case, one 
can discard the time component and view Mc as a compact subset of TM. Then Mather's 
theorem implies that is a Lipschitz graph for the projection vr : TM — ^ M. To see this, 
suppose that two curves x{t) and y{t) in 7r(Mc) have a:(0) = y{s) for some s. Mather's 
theorem rules out immediately the possibility that s is an integer, unless x = y is a periodic 
orbit. For a general s, consider the curve z{t) = y{t — s). Then, z{t) = y{t — s) and, by 
time-invariance of the Lagrangian, {z{t),z{t)) is a solution of the E-L flow. It has same 
average action and rotation vector as (y, y) and hence it is also in Mc- But then z{0) = x{0) 
is impossible, by Mather's theorem, unless i(0) = y{s) = x{0). 

In the realm of twist maps, one can also deduce from Mather's theory the existence of 
many invariant sets that are graphs over the base and are made of minimizers (see also 
[21]). 

4.2. EXAMPLES 

One would hope that /? is, in general, strictly convex, i.e. each point on graph{f3) is an 
extreme point. This is true when M = S^, and Mather shows in [28] how his Lipschitz 
Graph Theorem implies the classical Aubry-Mather Theorem, by taking a E-L flow that 
suspends the twist map. The fact that Mc is a graph nicely translates into the fact that 
orbits in an Aubry -Mather set are well ordered. 

The graph of (3 is also strictly convex when L is a Riemannian metric on T^. This was 
known by Hedlund [19] in the 30's, albeit in a different language. It will also be an easy 
consequence of Theorem 5. In Section 6, we will show that this may not be true if one 
adds a potential term to the metric in the Lagrangian. 

Here we briefly describe three other counter-examples to the strict convexity of /?. 

Example 3 (Mane [25]) Take L : TT^ R, given by L{x,x) = \\x - X|p where X is 
a vector field on T^. The integral curves x oi X are automatically E-L minimizers since 
L = on these curves. Mane chooses the vector field X to be a (constant) vector field 
of irrational slope multiplied by a carefully chosen function on the torus which is zero at 
exactly one point q. The integral flow of X has the rest point q{t) = q, and all the other 
solutions are dense on the torus. The flow of X (and its lift to TT^ by the differential) 
has exactly two ergodic measures: one is the Dirac measure supported on {q, 0), with zero 
rotation vector, the other is equivalent to the Lebesgue measure on and has nonzero 
rotation vector, say lj (see [18], for more details on this flow). Mane checks that /3~^(0) 
(trivially always an Xc) is the interval {Xoo \ A G [0, 1]}, and that no ergodic measure has 
a rotation vectors strictly inside this interval. Thus the Mather set Mq is the union of the 
supports of the two measures. 

In the context considered here, this example is a little unsatisfactory because it is not a 
mechanical Lagrangian. We will give an example in Section 6 of a mechanical Lagrangian 
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on which displays a similar phenomena. The next two examples involve metrics, on the 
three-torus and on the surface of genus two, respectively. 

Example 4 (Hedlund-Bangert) Consider in the three nonintersecting lines given 
by the x-axis, the y-axis translated by (0, 0, 1/2) and the z-axis translated by (1/2, 1/2, 0). 
Construct a Z^-lattice of nonintersecting axes by translating each one of these by all 
integer vectors. Take a metric in M.^ which is the Euclidean metric everywhere except 
in small, nonintersecting tubes around each of the axes in the lattice. In these tubes, 
multiply the Euclidean metric by a function A which is 1 on the boundary and attains its 
(arbitrarily small) minimum along the points in the center of the tubes, i.e. at the axes 
of the lattice. Because the construction is periodic, this metric induces a Riemannian 
metric on T''. One can show ([6]), if A is taken sufficiently small, that a minimal geodesic 
(which is a E-L minimizer in our context) can make at most three jumps between tubes. 
In particular, a recurrent E-L minimizer has to be one of the three disjoint periodic orbits 
which are the projection of the axes of the lattice. Thus there are only three rotation 
directions that minimizers can take in this example, or six if one counts positive and 
negative orientations. In terms of Mather's theory, the level sets of the function P are 
octahedrons with vertices (±a, 0, 0), (0, ±a, 0), (0, 0, ±a) (we assume here that A is the 
same around each of the tubes). Since we arc in the case of a metric, one can check that 
P is quadratic when restricted to a line through the origin (a minimizer of rotation vector 
au! is a reparametrization of a minimizer of rotation a;). Hence a set Sc is either a face, an 
edge or a vertex of some level set {/? = 6}, and the corresponding is, respectively, the 
union of three, two (parameterized at same speed) or one of the minimal periodic orbits 
one gets by projecting the disjoint axes. Note that, instead of the function /3 of Mather, 
Bangert uses the stable norm. Mather's function /3 is a generalization of that norm. 

It is important to note that the nonexistence of minimizers of a certain rotation vector 
cj does not mean that there are no orbits of the E-L flow that have rotation vector to. 
For example, Mark Levi has shown the existence of orbits of all rotation vectors in the 
Hedlund example (personal communication). In addition, in our torus example (Section 
6), it is easy to see that there are (nonminimizing) orbits with rotation vectors in the 
excluded interval. The next example was brought to our attention by A. Fathi. 

Example 5 Take the metric of constant negative curvature on the surface of genus 2 (the 
two-holed torus) which has a thin neck between the two holes (see Figure 1). In this case, 
the notion of minimizing is just that of least length using the hyperbolic metric. With 
a and b as shown, the minimal measure for the homology class a + b will be a linear 
combination of the ergodic measures supported on Fq and F;,, where F^ and F^ are the 
closed geodesies in the homotopy classes of a and b, respectively. This is because any closed 
curve that crosses the neck will be longer than the sum of the lengths of Fa and F;,. Hence 
(a -|- 6, P{a + b)) cannot be an extreme point of graph{(5). 

This example illustrates the remarks of the introduction, namely, on hyperbolic man- 
ifolds the notion of M minimizers is not the correct one if one wants to show that all 
the dynamics of the geodesic flow of a metric of constant negative curvature are preserved 
under global perturbation. The right notion of minimality comes from the observation that 
there is a curve of least length in the homotopy class of ab. This notion is generalized to 
an asymptotic invariant on hyperbolic manifolds that associates an ergodic measure of a 
given dynamical system with its "rotation measure" . The rotation measure is the ergodic 
measure of the hyperbolic geodesic flow whose dynamics best mirrors that of the given 
ergodic measure (see [9]). 



Dynamical Stability in Lagrangian Systems 



17 




5. Autonomous Lagrangian Systems on the 2- Torus 

We now give a concrete application of Mather's theory to autonomous Lagrangians on 
the two torus. If the Lagrangian is mechanical, then for each E > Vmax get a Jacobi 
metric. Hedlund showed ([19]) that for each Riemannian metric on the two torus there are 
minimizing geodesies in all directions. Applying this to the Jacobi metric we see that there 
will be minimizers in all directions on each energy level. However, the parameterization of 
these minimizers as solutions of the E-L flow is different than as Jacobi geodesies. Thus 
in applying Hedlund's result we must examine how the minimizers on the various levels 
fit together in terms of their rotation vectors. This can be done nicely in the framework 
of Mather's theory which also allows us to consider a wider class of Lagrangians. 

Recall that the sets Xc are projections on the rotation vector space of flat domains in 
graph{(5). 

Theorem 5 For any autonomous Lagrangian system on which satisfies Mather's hy- 
pothesis, a set Xc is either a finite interval in a line through the origin or a point. If Xc is 
an interval, either it contains 0, or it is supported by a line of rational slope. In addition: 

(i) There are minimal ergodic measures (and hence recurrent E-L minimizers) of all 
rotation directions. 

(a) If V has irrational direction, then there is a positive number K(v) such that \\v\\ > 

K{v) implies that there is a minimal ergodic measure (and hence a recurrent E-L 

minimizer) with rotation vector v. 
(Hi) In each rational direction, represented by a vector v, there is a sequence of positive 

numbers A„ — >■ oo and a sequence 7„ of periodic E-L minimizers such that 7„ has 

rotation vector XnV. 

(iv) If the Lagrangian is of the form L{x,x) = — V{x), where \\ ■ \\ is the Euclidean 
metric,^ then the origin is not in the interior of any Xc and the Xc which have the 
origin as endpoint have length less than 2y/Vmax — ymin- In particular, K(v) in (ii) 
can be taken to be 1\/Vm,ax — ^min- 

(v) The support of minimal measures can be either a point, a closed curve, the suspension 
of a Denjoy example (a lamination), or the whole torus. 

Proof By Mather's Lipschitz Graph Theorem and the remark after it, any set Mc projects 
injectively to T^. It thus yields a flow on compact set X C T"^. Most of the theorem is 
a consequence of the fact that there are severe restrictions on what X and its flow can 
be. We review some of these well known results. They follow easily from, for example, the 
techniques of [13] ( c/[34]). 

^Clearly, similar results can be obtained for other metrics 
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If X C is a compact set and there is a continuous flow (f)t dcflncd on X, then 
the support of a c/if-invariant ergodic measure on X can only be (1) a point, (2) a closed 
orbit, (3) the whole torus, or (4) homcomorphic to the suspension of the minimal set in a 
homeomorphism of the circle that is a Denjoy counter example (and it must be embedded 
in in the obvious way). 

In cases (1) and (2) there is a unique invariant probability measure on the set and so 
the rotation vector exists and is equal for all points. Case (1) happens only if the rotation 
vector is 0. Case (2) can happen only if the rotation vector is zero for homotopically trivial 
curves and a nonzero v G Q^, if the curve is homotopically nontrivial. Cases (3) and (4) 
are more subtle. The set X can contain a fixed point, but in any event the rotation vectors 
of ergodic measures supported on X can take on at most two values, zero and a vector 
with irrational slope. 

Of particular importance in what follows are the following three consequences. First, 
there cannot be two invariant measures in X whose rotation vectors are in different direc- 
tions, i.e. 7^ Second, if an ergodic measure has a nonzero rotation vector v G Q^, 
then it is supported on a closed orbit. Finally, X cannot support two ergodic measures 
with nonzero rotation vectors that have the same irrational slope. 

To prove the first paragraph of the theorem, recall that by Mather's theory, the end- 
points of each X^ corresponds to a minimal ergodic measure with that rotation vector. 
Next observe that any Xc is a finite interval on a line through the origin or a point, be- 
cause as noted in the previous paragraph, no Mc can contain orbits with different nonzero 
rotation directions and each Mc can contain ergodic measures that yield at most two ro- 
tation vectors. (The fact that Xc is finite also follows from the superlinearity of (3.) If we 
fix an ui with irrational slope, then if both endpoints of the interval X(. are nonzero, they 
correspond to at least two distinct ergodic measures with two irrational rotation vectors 
of same direction but different lengths, another impossibility. 

To prove (i), take any G R^. Using what we have just proved, any set Xc to which 
u) belongs is a finite interval with endpoints in the same direction as uj. These endpoints 
are the rotation vectors of two minimal ergodic measures. For (ii), note that if u has 
irrational direction, and if ||a;|| is large enough, we have just seen that Xc can only be a 
point. Hence there is a minimal ergodic measure of rotation vector u). As for (iii), note 
that the previous paragraph implies, in each rational rotation direction the existence of 
sequences A„ — > oo and minimal ergodic measures /x„ with p(/Ltn) = XnV. As noted at the 
beginning of the proof, such an ergodic measure necessarily comes from a periodic orbit. 

We now prove part (iv). We first show that /?(0) is an isolated global minimum, thus 
proving that cannot be in the interior of any Xc- In the case of autonomous mechanical 
systems, it is easy to see that a measure that gives the absolute minimum for A is the 
Dirac measure on a fixed point x at which V{x) = Vmax- This measure has rotation vector 
and action Vmax- 

Recall, moreover, that from Formula (2) in Section 4, we have the estimate f5{uj) > 
^\\uj\\'^ — Vmax for such systems. This implies that the graph of /3 above an Xc that contains 
zero cannot be horizontal, and that (0, Vmax) must be an isolated minimum of f3. 

Finally, we prove our estimate for the length of XcS which have as one of their 
endpoints. On an energy level E > Vmax the Jacobi metric \/E — V)ds is a Riemannian 
metric on T^. It is not hard to see that the Jacobi minimizers and the E-L minimizers 
for Le = ^{E — V)\\x\\'^ form the same set of curves on the torus (See [31], page 70). Le 
obviously satisfies Mather's hypothesis and so, by Part (i), there are Jacobi minimizers 
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of all rotation directions (we could also have invoked Hedlund [19], who proved this in 
the 30's). By Proposition 1, these Jacobi minimizers are also E-L minimizers for L = 
— V{x). In terms of minimal crgodic measures, after a time change we have found 
minimal crgodic measures of all rotation direction in each energy level E > Vmax- 
Let jj, be such a measure. Using Lemma 1, we have: 

^ {p{n)f - Vmax < J Ldn < J 2Vmindn = E- 2Vmin 

so pip) < \/2{E — 2Vmin + y-max)- This is true for the measures obtained from energy 
levels E arbitrarily close to Vmax- Thus we have shown that |jci;|| < 2y/Vmax — Vmin- 

Statement [v) is an immediate consequence of the facts at the beginning of the proof. 

□ 



6. Speed Defect: an Example in the 2-Torus 
6.1. THE MAIN EXAMPLE 

In this section, we exhibit an autonomous Lagrangian system on the torus which fails 
to have recurrent minimizers of all rotation vectors in a prescribed rational direction. This 
is a different kind of counterexample to a generalization of the Aubry-Mather theorem 
than Hedlund's example (one can also derive counterexamples in the realm of symplectic 

twist maps using examples in [2]). In terms of Mather's theory, we find in our system an 
X(. in a rational direction which contains an interval, and hence the function (5 of Mather 
is not strictly convex in this case, confirming the worst predictions of Theorem 5. 

Note that the Lagrangian in such an example could not be a Riemannian metric, since 
it is easy to show that Mather's function [3 for such systems is quadratic along lines through 
the origin of K^. Hence Theorem 5 above implies that, in the case of metrics, (3 is strictly 
convex. 

We first outline the features of our example. We will find, in a given energy level {h = 0), 
exactly two closed minimizers 71, 72 in the homology class (0, 1) G = ii'i(T^; Z) but with 
different rotation vector (0, pi), (0, p2)-, i.e. different average speeds. We then show that any 
periodic minimizer in that same homology class but in a different energy level goes strictly 
faster for higher energy, and strictly slower for lower energy than these two minimizers, 
leaving the segment between (pi,0) and {p2,^) empty of recurrent minimizers. 

On the torus = K^/27rZ^, consider a Lagrangian of the form: 



L(x,i) = ^^^\\xf-V{xi) (3) 
2 + cos 2x1 

Vixi) = — : (4) 

^ ^ 2 + sinxi ^ ' 

a{xi) = 2 + sinxi, (5) 

where ||i:|| denotes the usual, Euclidean norm on T^. The corresponding Hamiltonian is 
H{x,p) = ^ + V{x,). 

If there is a minimal ergodic measure in the homology direction (0, 1), then as noted 
at the beginning of the proof of Theorem 5, it must be supported on a periodic orbit. 
Furthermore, we claim that if 7 is a closed curve which is the projection of the support 
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of an ergodic, minimal measure in the (0, 1) direction, it has to be of the form xi = 
constant. Indeed, since H(x,p) is independent of X2, the Hamiltonian flow is invariant 
under translations in the X2 direction. Hence 7 is actually part of a one parameter family 
of periodic orbits that are translates of one another in the X2 direction. Obviously, all the 
orbits in this family have same rotation vector and action. Since 7 is the support of a 
minimal measure, all its X2 translates are as well, and they all belong to the same Xc- 
Because 7 is in the class (0, 1), it must have (at least 2) points at which = 0. If it is 
not so at all points, i.e. if xi is not constant along 7, it is easy to see that 7 will intersect 
transversally its X2 translates, contradicting (by Mather's theorem) the fact that these 
curves belong to the projection of the same X^. 

Our task is now to find all the possible periodic orbits of the Hamiltonian flow that 
project to curves xi = constant and to compute their rotation vector and action. Such a 
curve occurs when = xi = and hence (by setting pi(0) = 0), whenever 

. / ^ dH cx' ,, ,.9 ^ 

Simplifying this equation by fixing an energy level H = E and using = 2a{E — F), 
one gets: 

= cosa;i(2 + sina;i)(4sinxi — E) 

which has the solutions: 



Xl = |, (6) 

Xl = y, (7) 
4sinxi = E, (8) 

with some restrictions on what the range of E is in each case, which we will deal with 
later. 

We now find the curves that the corresponding ergodic measures trace in the (||a;|| = 
/>, A)- plane, and show that the third curve is always in the union of the epigraphs of the 
two first ones, and hence cannot be a minimal measure (the epigraph of a function is the 
set of points above its graph) . 

In the energy level {H = E}, the Lagrangian is given hj L = E — 2V. Since it is 
constant along the curves xi = constant, the average action of the corresponding ergodic 
(probability) measure n is 

A(/x) = J Ldfi = L J dfi = L 

From L = E — 2V, we can also get the rotation vector (0, p) of an xi = constant curve: 

. ...r a(xi) 2 X ^ 2 2(E-V(xi)) 

L + 2V = -^\xf + V(xi) = E^ \xf = — , , 

2 a{xi) 



In particular, the speed at which a curve xi = constant is traversed by the E-L flow is 
constant. The rotation vector of such a curve is: 
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a 



We found it more convenient to first look at the curves (c, A) formed by these measures, 
where 

a 

We denote by {ci{E) , Ai{E)) the curve corresponding to the i^^ family of solutions, i = 
1,2,3. 

• If xi = ^, Ai{E) = L = E - 2F(f ) =E + l, whereas ci{E) = ^ Jif^ = 1{E + i). 
Hence: 

Ai(c) = 3c+^. 

One can verify that the curve xi = | is indeed a solution for E in [—1/3, oo), which gives 
c G [0,oo). 

• If xi = ^, ^2(£^) = + 2, C2(-E) = E + 1 and hence 

-42(c) = c + 1. 

One can verify that xi = ^ is a solution for in [— 1, oo), which corresponds to c G [0, oo). 

• If £■ = 4sina;i, we get, by replacing sinxi = E/4 and cos 2x1 = 1 — 2sin^xi = 
l-2{E/4f: 

^E^ + 16E + 24: 
C3{E) = 2- 



AsiE) 



. (-g + 3) 
'(^ + 8)' 



We will be content to say that the curves £' = 4 sin xi correspond to solutions for E in 
some interval included in [—2,4] (one can check that — 2 < Vmin < ~1)) which we do not 
need to find, as we will see. 

Instead of trying to find ^3(0) exphcitly, we compute A3{E) - (csiE) + 1) = ^^(^+8)^^° 
which is positive for all E, thus showing that the curve {c,As{c)) is above (c, ^i(c)). The 
relation is trivially the same for the corresponding (p, A{p)) curves and hence the measures 
corresponding to the curves E = 4sinxi can never be minimal. 

Replacing c by 27r^/3^, we get: 

Ai{p) = 67rV + ^ 
A2{p) = 27rV + l, 

the graphs of which are two parabolas crossing at p = with Ai(p) < A2{p) for 

p G [0, A2 < Ai after that. Note that the union of the epigraphs of these curves 
is not convex. Remember that the graph of Mather's function /3 is the boundary of the 
convex hull of the points (p(p),A(p)) for all possible ergodic measures p . To form the 
lower boundary of the convex hull for the epigraphs of (p, ^^(p)), i = 1,2 (which, as we 
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O.'os ' ' o!l ' ' 0.'l5 ' ' 0^2 ' ' .25 ' ' o'.3 ' ' 0.35 



Figure 2. The graph of /3 in the (0, 1) direction 

showed, necessarily contains the epigraph of {p,A^(p)), we look for the line tangent to 
both these graphs. It is easy to check that it joins the point Qi = (^^^) 2/3) of the graph 

of Ai to the point Q2 = 2) of the graph of A2. Hence, any measure corresponding to 

points on the graphs of Ai or A2 above the segment of line between Qi and Q2 fails to be 
minimal. Instead, if Hi are the measures corresponding to Qi, i = 1,2, the minimal measure 
with rotation vector in the interval I = ]- , —i-l will be a convex combination of ui 

and /i2- The set for any oo = (0, p) with p G / is thus an interval, and the corresponding 
Mather set Mc is the union of the support of /ii and fi2, i-G. 2 loops, with rotation vectors 
(0, ) and (0, '^~)- Since any recurrent E-L minimizer with rotation number u> belongs 
to the Mather set Mc, it follows that there can't be any recurrent minimizer of rotation 
vector (0, p) with p in the interior of /. This is the gap we announced. 

Remark 4 It is easy to check that both pi,p2 have their support in {H = 0}, and 
that {p,Ai{p)) for p G [0, ^^^ ] corresponds to measures supported in £^ < whereas 

(p, A2{p)),p> corresponds to £^ > 0. 

It is instructive to think of the example in terms of the Jacobi metric. For E = 0, this 
metric is \/2cos(2xi)||x||, which has exactly two minimizers (with respect to the metric) 
in the homology class of (1, 0), namely, xi = 7r/2 and xi = 3tt/2. Of necessity these closed 
curves have the same Jacobi length, but when they are parameterized as solutions to the 
E-L equations, they have different speeds and thus different rotation vectors. 

As the energy E is swept through there is a "minimality exchange". Namely, when 
£' < 0, a;i = 7r/2 is the Jacobi minimizer and corresponds to a minimal measure for the 
E-L flow. When > 0, xi = 3it/2 has these properties. At E = they share these 
properties, but as noted above, they have different rotation vectors for the E-L flow. 

6.2. A E-L MINIMIZER WHICH IS NOT A JACOBI MINIMIZER 

Using the example above we construct a E-L minimizer that is not a Jacobi minimizer on 
its energy level, thus showing that the converse of Proposition 1 in Section 2.4 is false. 
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Let ctE he the parameterization of the curve xi = it/ 2 with the property that daE 
is a trajectory of the E-L flow with H = E e [1/3, oo). Similarly, /3e corresponds to the 
curve = 37r/2 for £^ G [— 1, oo). We shall use a and (3 to refer to these curves without a 
specific parameterization. 

If He and r]E are the ergodic invariant probability measures supported on daE and d^E, 
respectively, then we have that for E' < and E > 0, respectively, he and r)E are minimal 
measures. Further, /ig and rjQ arc the minimal measures corresponding to the rotation 
vectors pi = (0, and p2 = (0, -^) and for any p = (0,p) with < p < 

the minimal measure is a linear combination of pQ and rjQ. In particular, if = and 

P3 = (0,^3) = , then the minimal measure with rotation vector ps is ^^2±Z2o_ 

Now for each N eN, let 77V be the E-L minimizer satisfying ^j^{—N/2) = (tt, —Np^/2) 
and ^]\f{N/2) = {tt, Np3/2), and so p{j'n) = ps- If wc let be the probability measure 
uniformly distributed with respect to time on d'jj\f, then the results in section 2 of [27] 
says that as ^ 00, will converge weakly to the minimal measure with rotation 
vector pz and so ctat —>■ ^^2±Z2o jj^ particular, by choosing N large, we can make H{'jn, Jn) 
arbitrarily close to zero, and we can insure that 7jv spends about half its time near ao 
and the other half near Pq. For each N, let ijsi be the path that is the union of 

[(7r/2, (tt, ^)], [(7r/2, ^), (7r/2, )], and [(7r/2, ), (tt, )], 

where we use square brackets to indicate the straight segment connecting two points in 
the plane. Thus ^jv goes directly left from the top endpoint of 7iv, then travels down a 
and then directly across to the bottom endpoint of 7jv. 

Now we assume that for all A/", 77V is a Jacobi minimizer for its energy £Jjv and obtain 
a contradiction. For N large, we have that |£^jv| is small. Now if £^jv < 0, since is the 
minimizer for those E, it is clear that ^at has lesser Jacobi length than 77V for energy Ef^ 
because 77V must spend half of its time near /?. By considering a curve analogous to 
but going to the right and then down /3, we find that E^ > is impossible also. 

Remark 5 As noted in Section 3.3 and Example 4 in Section 4.2, when one attempts 
to take the limit of minimizers 77V without adequate control, it can happen that the 7jv 
converge to something that is not in the same direction as each ^fq. The example of this 
section illustrates another mechanism by which the limit argument can fail. It can happen 
that the 7Ar converge to something in the correct direction and the speed of the limiting 
measures average to the correct speed, but the various ergodic components of the limit 
measure have either greater or lesser speed than the desired one. This gives rise to the 
speed gap one has to allow in Theorem 1 and Theorem 5 . 
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Note Added in Proof 

Just prior to the publication of this paper we became aware of the elegant paper "On 
minimizing measures of the action of autonomous Lagrangians" , by M. J. Dias Carneiro 
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{Nonlinearity, 8, 1077-1085, 1995) which contains much that is relevant to Sections 5 
and 6 of this paper. In particular, Dias Carneiro proves that in the autonomous case, 
the function /3 (defined in Section 4.1) has a derivative in radial directions, calculates the 
derivative, and as a corollary, shows that Mather sets are always supported on a single 
energy level. These results simplify the required calculations in the example of Section 6. 
In addition, Dias Carneiro makes a remark that contains part of Theorem 5 and gives a 
result connecting the minimal measures of the Jacobi metric with those of the E-L flow. 

In addition, R. Iturriaga showed us notes of Mane which contain 11 beautiful theorems 
about autonomous Lagrangian systems. His students are providing the proofs of these 
results. 

We dedicate this paper to the memory of Ricardo Maiie. 
References 

1. Abraham, R. and Marsden, J.E. (1985) Foundations of Mechanics, Addison- Wesley 

2. Arnaud, M.C (1989) Sur les points fixes des diffeomorphismes exacts symplectiques de T" x M", C.R. 
Acad. Sci, Paris 309, pp. 191-194 

3. Arnold, V.I. (1980) Mathematical methods of Classical mechanics, Springer- Verlag 

4. Aubry, S. and Lc Daoron, P.Y. (1983) The discrete Frenkel-Kontorova model and its extensions I: 
Exact results for the ground states, Physica 8D, pp. 381-422 

5. Bangort, V. (1988) Mather sets for twist maps and geodesies on tori, Dynamics Reported 1 (eds. 
Kirchgrabcr, U. and Walther, H.O.) John Wiley 

6. Bangcrt, V. (1989) Minimal geodesies, Ergod. Th. & Dynam. Sys., 10, pp. 263-286 

7. Bernstein, D. & Katok, A. (1987) Birkhoff periodic orbits for small perturbations of completely inte- 
grable Hamiltoniam systems. Invent. Math. 88, pp. 225-241 

8. Boyland, P. (1994) Topological methods in surface dynamics Topology and its Applications 58, pp. 
223-298 

9. Boyland, P, Asymptotic dynamical invariants on hyperbolic manifolds, in preparation 

10. Boyland, P. and Golc, C. (1995) Lagrangian systems on hyperbolic manifolds, preprint 

11. Bialy, M. and Poltorovitch, L. (1992) Hamiltonian systems, Lagrangian tori, and Birkhoff's theorem 
Math. Ann. 292, pp. 619-627 

12. Denzler, J. (1987) Mather sets for plane Hamiltonian systems. J. Appl. Math. Phys. (ZAMP) 38, pp. 
791-812 

13. Prank, J. and Misiurewicz, M. (1990) Rotation sets for toral flows, Proc. A. M.S. 109, pp. 243-249 

14. Puller, P. B. (1965) On the surface of section and periodic trajectories. Am. J. Math., 87, pp. 473-480 

15. Ghys, E. and de la Harpe, P. (eds) (1990) Sur les groupes hyperboliques d'apres Mikhael Gromov, 
Birkhausor 

16. Golc, C. (1994) Periodic orbits for Hamiltonians in cotangent bundles Trans. A. M.S. 343, Number 
1, pp. 327-347 

17. Gromov, M. (1977) Three remarks on geodesic dynamics and fundamental group, preprint, SUNY at 
Stony Brook 

18. Handel, M., Periodic point free homeomorphisms of T^, Proc. A. M.S. 107, pp. 511-516. 

19. Hedlund, G.A. (1932) Geodesies on a two-dimensional Riemannian manifold with periodic coefficients. 
Ann. of Math 33, pp. 719-739 

20. Herman, M. (1989) Inegalites "a priori" pour des tores Lagrangiens invariants par des 
diffeomorphismes symplectiques, Publ. Math., I.H.E.S. 70, pp. 47-101 

21. Katok, A. (1992) Minimal orbits for small perturbations of completely integrable Hamiltonian sys- 
tems, Twist mappings and their applications, I.M.A. 44 (McGehee, R. and Meyer, K. eds) Springer- 
Verlag 

22. Klingenberg, W. (1982) Riernannian Geometry Walter de Gruyter 

23. MacKay, R.S., Meiss, J.D. and Stark, J. (1989) Converse KAM theory for symplectic twist maps, 
Nonhnearity 2, pp. 555-570 

24. Mane, R. Ergodic theory and differentiate dynamics. Springer- Verlag 

25. Mafie, R. (1991) Global variational methods in conservative dynamics, 18° Coloquio Brasileiro de 
Matematica, IMPA 

26. Mather, J. (1991) Variational construction of orbits of twist difFeomorphisms, ,Jour. of A. M.S. 4, pp. 
207-263 

27. Mather, J. (1991) Action minimizing invariant measures for positive definite Lagrangian systems. 



Dynamical Stability in Lagrangian Systems 



25 



Math. Z. 207, pp. 169-207 

28. Mather, J. (1993) Variational constructions of connecting orbits, Ann. Inst. Fourier, Grenoble 43, 5, 
pp. 1349-1386 

29. Meiss, J. D. (1992) Symplectic maps, variational principles, and transport. Rev. Mod. Phys. 64, No. 
3, pp. 795-848 

30. Denvir, J. and MacKay, R.S. (1995) Consequences of contractible geodesies on surfaces, preprint 

31. Milnor, J. (1969) Morse theory, Princeton University Press 

32. Morse, M. (1924) A fundamental class of geodesies on any closed surface of genus greater than one, 
Trans. Math. Soc. 26, pp. 25-60 

33. Moser, ,J. (1986) Monotone twist mappings and the calculus of variations. Erg. Th. & Dynam. Sys. 
8*, pp. 199-214 

34. Walsh, J. (1995) Rotation vectors for toral maps and flows: a tutorial. Int. Joum. of Bif. and Chaos 
5, pp. 321-348 



